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This is the first of two papers dealing with single server systems. Two 
subjects are discussed in the present paper, (i) relations between such items 
as the probability of loss, probability of no delay, and the average number of 
customers served in a busy period, and (ii) the statistical behavior of a 
single server system in which no waiting for service is allowed. 

I. INTRODUCTION 

A typical but rather homely example of the systems considered here is 
a barber shop in which there is only one barber, the "single server." 
Let N — 1 be the number of chairs provided for customers waiting for 
service so that the "capacity" of the system is N. When the shop is 
full, one customer is being served and N — 1 are waiting. A prospective 
customer (a "demand" for service) arriving when the shop is full is 
turned away and is said to be "lost." If the shop is not full, he waits and 
is eventually served. 

The demands arrive at an average rate of a per unit time. The server 
would serve b customers per unit time (on the average) if he were to work 
steadily. It follows that the average interval between arrivals is 1/a and 
the average service time is 1/6. It is assumed that the rates a and b do 
not change with time. 

The first part of the paper is concerned with several quantities of 
interest, including the fraction L of demands lost and the average length 
of the busy periods, i.e., the periods during which the server is continu- 
ously busy. The values of these quantities are expressed in terms of a, 
b, and two other quantities p and r. Here p is the probability that the 
server is idle (at an instant selected at random) and t is the average 
duration of an idle period. Both p and r depend upon N and upon the 
probability laws governing the arrivals and service times. However, only 
the simplest cases of this dependence are mentioned in the first part, of 
the paper. The results are summarized in Table I. 
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The second part of the paper is concerned with the single server "loss 
system" in which no waiting is allowed (N = 1). The input is now as- 
sumed to be "recurrent," i.e., the distances between arrivals are inde- 
pendent and have the general distribution function A(t). The service 
times have the distribution function B(t) and are independent of each 
other and of the arrivals. The functions A(t) and B(t) are such that the 
interarrival distances and service lengths have the respective average 
values 

aT 1 = f [1 - A(t)]dt 

(1) 
b' 1 = [ [1 - B(t)] dt. 

Further conditions are imposed on A(t) and B(t) in the course of the 
derivations. The principal items of interest are (i) the loss L, (ii) the 
probabilities po and pi that the server is idle or busy, respectively, at a 
time selected at random, and (Hi) the distribution of the lengths of the 
idle periods. The expression for the loss is equivalent to one obtained in 
a different manner by F. Pollaczek. 1 A discussion of how the loss in- 
creases with the variability of arrival has been given by P. M. Morse. 2 

Results for a Type I counter, i.e., one in which the registration of an 
arrival is followed by a "dead time," may be applied to the single server 
loss system by identifying the dead time (supposed variable) with the 
intervals the server is busy. For example, results obtained here are 
closely related to some for counters given by R. Pyke" 1 and earlier workers 
to whom he makes reference. 

The infinite capacity system (N = °o ) is discussed in a companion 
paper. 4 In this case there is no loss, and attention is focused on the 
distributions of the waiting times and busy period lengths. Some of the 
results of Section II of the present paper find application there. 

I am indebted to John Riordan for many helpful discussions on the 
subject matter of these two papers and for numerous improvements in 
presentation. 

II. AVERAGES OBTAINED FROM FIRST PRINCIPLES, GENERAL INPUT AND 
LIMITED CAPACITY 

Several important averages associated with a single server system may 
be readily obtained by considering its behavior over a long period of 
time T, statistical equilibrium being assumed, and then letting T — > =» . 
The input and service are assumed to be general with respective arrival 
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and service rates a and b. The capacity of the system is N. so that de- 
mands arriving when JV — 1 are waiting are lost. 

The probability po that the server is idle (at an instant selected at 
random) and the average duration t of an idle period are supposed given. 
For Poisson input, i.e., one in which the probability that a demand will 
arrive in the interval t, t + dt is adt + 0(dt ) (irrespective of the arrival 
times of the other demands), the value of t is 1/a. 

Let v(T) be the number of arrivals in T and let e be a given, arbi- 
trarily small, positive number. The input and service are assumed to be 
such that the probability of 1 — e < v{T)/aT < 1 + e may be made 
as close to unity as desired by choosing T large enough. For recurrent 
input this restriction is satisfied, by virtue of the law of large numbers, 
when the first integral in ( 1 ) gives a finite value for a -1 . We shall refer 
to the above inequality by saying that the number of arrivals in the 
long interval T is "equal" to aT. Similar statements made below regard- 
ing total idle and busy time, number served, etc., in the interval T are 
to be interpreted in a similar way. 

The total idle time of the server is Tpa , the total busy time is T — Tp , 
and the total number served is b(T — Tp ). The number lost is the 
number of arrivals less the number served, and the fraction of arrivals 
lost is 

L = [aT - b(T - Tp )]/aT = 1 - (1 - p )p _1 (2) 

Avhere p = a/b. L is the "probability of loss." Loss occurs only when the 
waiting room is filled, i.e., when the system is in state N (probability 
Pn and average duration t n ). The number of times state N occurs is 
PsT/t n , and the average number of demands lost during such a period 
is r N [a — 6(1 — Pg)]/pn . For exponential service the service time 
lengths have the distribution function B(t) = 1 — e~ bl , and the value of 
tjv is 1/6. 

The number served in T without delay is equal, to within one, to the 
number of idle periods. The number of idle periods is poT/r; hence po/ar 
is the probability W(0) that a demand will be served upon arrival. Thus, 

W(0) = po/ar. (3) 

Because a demand is either served without delay, delayed, or lost, the 
probability of delay is 

1 — L — ( po/ar) = (1 — po)p _1 — (po/ar). 

Since busy and idle periods alternate, the number of busy periods is, 
to within one, the number of idle periods, namely poT/r. Dividing this 
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into the total busy time T — Tp gives t(1 — po)/po for the average 
length of a busy period. Similarly, the total number served gives br( 1 — 
po)/po for the average number served in a busy period. 

When N is infinite, statistical equilibrium requires p < 1, and' all 
demands are eventually served. Equating the number of arrivals a? 1 to 
the number served b(T — Tpo) gives p = 1 — p- 

These results and the forms they assume for N infinite and t = 1/a 
(Poisson input) are summarized in Table I. The results for Poisson 
input are well known. 

When the capacity of the waiting room is infinite, the average number 
n of demands served in a busy period is 



n = l/W(0). 



(4) 



In this case L is zero, demands are either delayed or not delayed, and 
W(0) becomes the probability of no delay. This follows from (3), Table 
I, and the fact that both sides of (4) are equal to ar/(l — p). For limited 
capacity 

n = (1 - L)/W(0) (5) 

which follows from Table I. Here W(Q) is the chance that a demand 
chosen at random will be served upon arrival. 

III. SINGLE SERVER LOSS SYSTEM 

In a single server loss system any demand arriving when the server is 
busy is lost. The system capacity N is 1, there is no waiting line, and 
every busy period consists of a single service interval. As mentioned in 
the introduction, the input is assumed to be recurrent, the service is 



Table I — Stationary Averages for a General System 





Limited Capacity 
General Input 


Infinite Capacity, po = 1 — p 


Average 


General Input 


Poisson Input. 

T = a"' 


Rel. number lost 
Rel. number not de- 
layed 
Rel. number delayed 
Length of busy period 
Number served in 
busy period 


1 - (1 - po)p- 1 
Po(fflr) ' 

(1 — po)p~' - po(aT-) -1 
t(1 - p )po 1 
6t(1 — po)po l 




(1 - p)(ar)-' 

1 - (1 - p)^)- 1 
r P (l - p)- 1 
ar(l - P )-i 




1 - p 

(b - a)- 1 
(1 - p)- 1 



Notation: a = arrival rate, b = service rate, p = a/b, t = average length of 
idle period, po = fraction of time server is idle. 
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general, and the principal items of interest are (i) the loss L, (ii) the 
probabilities p and pi that the server is idle or busy, respectively, at a 
time selected at random, and (Hi) the probability q(u) du that the 
length of an idle period will lie between u and u + du. 

3.1 Values of L, p () and p\ 

For general input and service, the results of Section II for N = 1 and 
the relation I — p = pi show that 

L = 1 — pip -1 , po = brpi , p = a/b (6) 

where the second relation is obtained by equating the average busy 
period length t(1 — po)/po to the average service length 1/6, t being 
the average idle period length. 

For Poisson input, t = 1/a and (6) gives 

Pi = PPo, p = (1 + p)~\ L = p x = p/(l + p). (7) 

In this case, the loss is independent of the service distribution, a property 
also possessed by the many-server loss system for Poisson input. 

For recurrent input and general service, the ratio L/(\ — L) of the 
number of demands lost to the number served is equal to the expected 
number of demands arriving while the server is busy serving one de- 
mand, i.e., during one service interval. Thus 

— ^ = fnTT)dB(t) (8) 

where n(t) is the expected number of arrivals during a service of length t 
(not counting the one starting the service) and B(t) is the service time 
distribution function. As t becomes large, n(t) is 0(at) and the integral 
converges because the average service time is finite. 
It will be shown that 

W) =±-.r" .f* )e ". d *. U>Q) (9) 

2in J c -,» [1 — a(s)]s 

where a(s) is the I.aplace-Stieltjes transform of the distribution function 
A(t) for the separation between arrivals: 

a (s) = f c- 8 'd,A(t). (10), 

It is assumed that a(s) and t are such that the integral in (9) converges. 
Equation (9) is but one of a number of similar results; see, for example, 
D. R. Cox and W. L. Smith. 6 
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To obtain (9), note that the service starts with an arrival, and the 
probability that n or more additional arrivals will occur in the ensuing 
interval of length t is the probability that 

S n = Z! + X 2 + ••• +Z„ ^ t. (11) 

Here X, is the separation between arrivals i — 1 and i. Since the X,-'s 
are independent and have the distribution function A(t), the rules for 
determining the distribution of the sum of n random variables may be 
applied to find the chance that S n ^ t. In particular, the Laplace trans- 
form of the probability density for <S„ is [a(s)] n , and the Laplace trans- 
form of Prob [S n g t] is [a(s)] n /s. 

The probability of exactly n arrivals in an interval of length t which 
starts just after an arrival is 

P„(t) = Prob[£„ ^ t) - Prob[& + i ^ t] 

J_ r+ ix . , .. . (12) 



= J-. \ C [1 - aismctisWs-V 1 ds. 



Multiplying P n (t) by n, noting that | a(s) ( < 1 on the path of integra- 
tion, and summing from n = to n = °o then gives (9). 

When n(t) is known as a function of t, either from (9) or otherwise, 
and is used in (8), the result is an equation which may be solved for the 
loss L. When L is known, (6) gives 

Pl = P (l-L), Pd-1-Pi. (13) 

A few examples follow. 

i. Poisson input. Here 1 — A(t) = e~ at and a(s) = a/ (a + s). Then 

w (0 = o - : / as~ 2 e st ds = at 

, L _ = / atdB(t) = a/b = P 
1 — L Jo 

and the results given in (7) are again obtained. 
it. Arrivals spaced 1/a apart. By inspection 

nit) = n, na~ l < t < (n + l)a -1 

and (8) becomes 

L 
1 - L 



?»/.r" <a( ' ) -?-[ B ( !L T i )- B G)]- (u) 



Hi. When a(s) is 0(l/s) as s — > c ± i°o , as it is when the probabihty 
density A'(t) = dA(t)/dt exists and is of bounded variation, and when 



(15) 
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B'(t) exists and is 0(e~"), e > 0, as t -> <», substitution of (9) in (8) 
gives 

1 - L h 2iri i c -.» [1 - a(s)]s 

= ir^a(sM-s)ds (0<c<e) 
27rt «/«-,•* [1 — a(s)\s 

where 

(3( s ) = I e - sl dB{t). (16) 

In (15) the singularities of (5{—s) lie to the right of the path of integra- 
tion and those of a(s)/s[l — a(s)] to the left. The conditions imposed 
on a(s) and B'(t) are sufficient to ensure the absolute convergence of 
the double integral and hence to justify the inversion of the order of 
integration. The result (15) is equivalent to one obtained by F. Pollac- 
zek 1 by a different method. 

iv. Exponential service. Here, 1 — B(t) = e~ ' and P(—s) = 
6/(6 — s). Substituting this value of /3( — s) in (15), closing the path of 
integration by an infinite semicircle on the right, and evaluating the 
residue at the pole s = b gives 

L " ,/,i L = a(b), Pl = p[l - *(&)]. 



1 - L 1 - «(&)' 

This expression for pi is a special case of the stationary state probabilities 
(determined by both F. Pollaczek 6 and L. Takacs 7 ) for the many server 
system with recurrent input and exponential service. 

v. Let a(s) satisfy the same condition as in example Hi and in addi- 
tion, suppose that (9) may be written as 

m = w) + ± r" ,? s)e ", d ° (oo) 

2irl J-c-i« [1 — a(s)]s 

where the only singularity of e 8t a(s)/s[\ — a(s)} to the right of Re(s) = 
— c(s finite) is a double pole at s = with residue R(t). Then (8) gives 

-Lj = r R( t) dB (t) + ' r" a -^fi^i . d7) 

1 — L Jo 2in J-c-ix> [1 — a(s)]s 

It should be noticed that not all cases can be handled by (15) and (17). 
An example to the contrary is furnished by 

A(t) = 1 - (1 +a0 -2 

(18) 
B(t) = 1 - (1 + 60 " 
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where the averages cT and 6 _1 exist but the corresponding variances are 
infinite. In this case 

a(s) = 1 — x + xV / e - "*/ -1 dy (x = s/a). 

As | s I -> oo in Re(s) ^ 0, a(s) is 0(s _1 ). Near s = 0, 

a(s) = 1 — saT 1 — s 2 a~ In s + O(s'). 

Since /3(s) is similar to a(s), both £( — s) and a(s)/s[l — a(s)] have 
branch points at s = 0. Equation (17) fails in this case because s = 
is not a double pole. Equation (15) fails because it is impossible to draw 
a path of integration which separates the singularities of 0(—s) from 
those of a(s)/s[l — a(s)]. When A(t)&nd B(t) are given by (18) it 
seems necessary to work directly with (8) and (9), or else use some sort 
of limit. 

3.2 Lengths of Idle Periods 

Now turn to the distribution of l b and U , the lengths of a busy period 
and the following idle period, recurrent input and general service being 
assumed. The probability that h ^ t is simply B(t), the service length 
distribution function. The distribution of the idle period length /, is more 
complicated. Its average length is, from (6), 

t = p /b Pl = a"'(l - L)~ x - /T 1 (19) 

where L is given by (8). The first step towards obtaining the probability 
q(u) du that u < U < u -\- du is to determine the conditional prob- 
ability q(u; t) du that u < U < u + du given l b = t. 

Consideration of the arrival patterns which give no arrivals in 
(t, t + u) followed by one in {t + u, t + u 4- du) leads to 

CO 

q(u; t) du = £ Pr [S n < t\ t + u < *S„+i < t + u + du] (20) 

71=0 

where So is and *S„ for n > is the sum (see (11)) of n interarrival 
intervals X,- . An expression for the joint probability density of S n and 
*S„ + i may be obtained by inverting its double Laplace transform 

ave exp [— rS„ — s»S n+ i] 

= ave exp [-(r + s)(X 1 + • ■ • + X n ) - sX n+1 ] 
= [a(r + s)] n a{s). 
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Integrating the density over the region < S„ < t, u + t < S„ +i < 
u + t + flu shows that the nth term, n > 0, in (20) is 

P^ 2 r iX dse° (u+l) a(s) r X [a(r + s)] n (e rl - Or" 1 dr (21) 

2m) 2 J C -i<>> ->c-ix 



(2«) 



where c > 0. 

Expression (21) holds only for n > 0. However, replacing the factor 
(e rt — 1) by e Tt gives an expression which holds for n ^ 0. Indeed, clos- 
ing the path of r-integration on the right shows that the integral of 
[a(r + s)]"r _1 is zero for n > 0. Closing it on the left shows that the 
integral of e rt r~ l is 2iri for t > and leads to the correct value for n = 0. 

Setting the modified form of expression (21) in the series (20) and 
performing the summation shows that q(u; t) du is equal to an expres- 
sion obtained by replacing [a(r + s)]"(e Tl — 1) in (21) by 
[1 — a(r + s)]~ 1 e rt . The joint probability density of h and U is 
q(u; t)B'(t) where B'(t) = dB(t)/dt. The probability density q(u) is 
the integral of q(u; t)B'(t) taken from t = to t — °o. Assuming B'(t) 
to be 0(e~ fl ) as / — * ■» and choosing the paths of integration c ± i°o 
so that < 2c < e makes the integral of B'(t) exp [t(s + r)] converge 
and have the value /3( — s — /•). Changing the variable of integration 
from r to z = r -f- s and, for convenience in writing (22), taking a(s) to 
be such that the path of integration for s may be shifted to — ij ± too 
(this implies that A'(t) is 0(e~'") as t — > » ) gives finally 

q(u) =[—~\f dse a(s) / r rTf7 r (22) 

\27iV J_,_,-w J,-t« [1 — a(z)](z — s) 

where u ^ and t\ is an arbitrarily small positive number. 

It may be shown that (22) reduces to a e~ au for Poisson input, as it 
should, and to 

q(u) = 6[1 - ot(b)]- 1 [ e- bv A'(u + v) do (23) 

for recurrent input and exponential service. Multiplying q{u) by 
exp (— s'u), integrating u from to °o, closing the path of integration 
for s on the right, and dropping the prime from s' shows that the Laplace 
transform of q(u) is 



-su -si, 

ave e = ave e 



i,r^p';»-fl t (24) 
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We also have 

-rl b -* 



avee 



... = i r* efr^) r«(«) - «wi * (25 ) 

2ire J,-,-* 2 — s [ 1- a(z) J 



which may be regarded as a double Laplace transform. In (24) and (25) 
the singularities of /3(— z) and (i(r — z) are supposed to lie to the right 
of the path of integration and the remaining singularities of the inte- 
grands to the left. 
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